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ABSTRACT 


/ 

For  many  multivariate  hypotheses,  under  the  normality  assump- 
ions,  the  likelihood  ratio  tests  are  optimal  in  the  sense  of 
aving  maximal  exact  slopes.  The  exact  distributions  needed  for 
implementing  these  tests  are  complex  and  their  tabulation  is 
limited  in  scope  and  accessibility .  In  this  paper,  a  method  of 
constructing  normal  approximations  to  these  distributions  is 
described,  and  illustrated  using  the  problems  of  testing  sphericity 
and  independence  between  two  sets  of  variates.  The  normal  approxi¬ 
mations  are  compared  with  well  known  competing  approximations  and 
are  seen  to  fare  well. 


Key  Words: 


Sphericity,  independence  between  two  sets  of  variates. 


<7,3 


♦Research  supported  in  part  by  the  Air  Force  Office  of  Scientific 
Research,  Air  Force  Systems  Command,  USAF  under  Grant  No.  AFOSR- 
77-3360.  7  The  United  States  Government  is  authorized  to  reproduce 
and  distribute  reprints  for  Governmental  purposes  notwithstanding 
any  copyright  notation  hereon. 


Accession  For 


f  7 


A 


*  — . -  - 


-1- 


1.  INTRODUCTION 

For  most  testing  of  hypothesis  problems  in  multivariate  anal¬ 
ysis,  under  the  normality  assumption,  several  reasonable  solutions 
of  comparable  merit  exist.  These  include  the  tests  resulting  from 
the  union-intersection  principle,  the  class  of  likelihood  ratio 
criteria  and  adhoc  statistics  such  as  Bar tlet t-Pillai  trace  for 
MANOVA.  The  Neyman-Pear son  theory  provides  some  information  on 
the  operating  characteristics  of  these  procedures,  but  does  not 
indicate  any  of  the  contenders  as  superior.  However,  as  demon¬ 
strated  by  Hsieh  (1979),  the  likelihood  ratio  tests  for  many  of  the 
multivariate  hypotheses  have  maximal  exact  slopes,  i.e.,  they  are 
asymptotically  optimal  according  to  Bahadur's  (1967)  method  of 
comparing  tests.  From  a  practical  standpoint  the  null  distri¬ 
butions  of  the  likelihood  ratio  statistics  or  of  their  competitors, 
are  of  crucial  importance.  These  distributions,  where  available, 
are  complex,  their  tables  are  generally  limited  In  scope  and  not 
often  accessible.  Moreover,  the  tabulations  concern  only  selected 
percentiles  and  are  inadequate  for  computing  the  p-values  needed 
in  practice.  The  pragmatic  approach  to  such  distribution  problems 
from  early  days  (e.g.,  Neyman  and  Pearson,  1931)  is  to  seek  reason¬ 
ably  accurate  and  convenient  approximations  to  the  distributions. 

The  principal  methods  of  approximating  a  likelihood  ratio  use 
the  fact  that,  in  large  samples,  its  distribution  is  approximately 
of  Pearson  type  I  form  and  that  of  its  negative  logarithm  is  of 
type  III,  i.e.,  chi-square,  form.  Nayer  (1936)  following  a  sug¬ 
gestion  by  Neyman  and  Pearson  (1931)  used  the  moments  to  approxi- 
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mate  the  percentiles  for  testing  the  homogeneity  of  variances  in 
this  manner.  Bishop  (1939),  on  the  other  hand,  obtained  empirical 
expressions  for  the  parameters  for  a  type  I  approximation  by  pass¬ 
ing  the  intermediate  stage  of  computing  the  moments.  Bartlett 
(1937)  pursuing  the  asymptotic  chi-square  character  of  a  negative 
multiple  of  loglikelihood  ratio,  pointed  out  by  Neyman  and  Pearson 
(1931) ,  used  moments  to  approximate  it  by  a  scaled  chi-square 
variable  for  samples  of  moderate  size.  This  approximation  deteri¬ 
orates  as  the  size  of  the  problem,  as  measured  by  the  dimension  of 
the  multivariate  normal  distribution  or  by  the  number  of  populations 
in  the  problem  increases,  or  when  the  effective  sample  size  is  small. 
A  comprehensive  investigation  of  various  approximations  was  con¬ 
ducted  by  Box  (1949),  in  which  he  introduced  new  widely  known  and 
used  asymtotic  chi-square  series  approximations  for  the  distribu¬ 
tions  of  likelihood  ratios.  Box  studied  his  series  approximations, 
in  the  context  of  two  multivariate  problems,  comparing  them  with 
the  exact  distributions  and  with  several  other  approximations 
including  one  based  on  the  F-dis tribution. 

The  purpose  of  this  essay  is  to  describe  a  method  for  con¬ 
structing  a  Gaussian  approximation  to  the  null  distribution  of  the 
likelihood  ratio,  and  to  demonstrate  its  efficacy  and  relevance  in 
testing  multivariate  hypotheses.  The  normal  approximation  is  out¬ 
lined  in  Section  2.  It  is  illustrated  using  two  common  multi¬ 
variate  problems,  namely  testing  independence  of  two  sets  of 
variates  and  testing  the  sphericity  hypotheses.  Section  3  con¬ 
tains  the  likelihood  ratio  statistics  for  the  two  problems  together 


with  current  approximations  for  their  null  distributions.  These 


approximations  are  then  numerically  compared  with  the  new  normal 
approximation  in  Section  4. 


2.  A  NORMAL  APPROXIMATION 
FOR  THE  LIKELIHOOD  RATIO  A 

Let  Y^»  .  ••»  Yq  be  a  sequence  of  asymptotically  normally 

distributed  nonnegative  random  variables.  The  convergence  of  the 

distribution  of  Yq  to  normality  can  be  accelerated  by  approximately 

symmetrizing  it  with  a  transformation  as  follows: 

Let  K  *  k  (n) »  r  »  1,  2,  .  ..»  denote  the  cumulants  of  Y  *  Y 
r  r  n 

and  suppose  that  and  k^/k^  m  $r  >  r  >_  2 »  are  bounded  as  n-*». 

Then  using  the  Taylor  series  it  is  easy  to  obtain  the  following 
asymptotic  expansion  for  the  expectation  E(Y/k^)^  of  a  power  of  Y 
as 


«>> .  i  +  +  tH 

K1  24tcJ 


+  3(h-3)$3}  +  0(k"3). 


(1) 


From  this  the  moment  of  (Y/ic^)h  can  be  obtained  by  substitut¬ 

ing  (rh)  for  h  in  (1).  The  following  central  moments  of  (Y/k^)*1 
are  then  obtained  in  a  routine  manner: 


P2(h) 
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P  3  (h)  -  ^{<^3  +  3(h-l)<fr*}  +  0  (ic~3)  , 

K1 
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(3) 


...  ... 


Since  Y  is  asymptotically  normally  distributed  as  n-^,  by  the  Mann- 
Wald  (1943)  theorem  so  is  an  appropriately  normalized  (Y/ic^)*1. 

This  convergence  to  normality  is  accelerated  if  h  is  chosen  so  that 
the  leading  term  in  the  expansion  (3)  for  y^(h)  vanishes.  This 
value  hQ  of  h  which  approximately  symmetrizes  (Y/tc^)**  is  obtained 
from  (3)  as 

hQ  -  1  -  IC1IC3/  (3K2>  * 
h 

The  distribution  of  (Y/k^)  may  be  approximated  by  the  normal  dis- 

f 

tribution  with  mean  p^(hQ)  an<*  variance  given  in  (1)  and  (2), 

respectively.  That  is, 

hn 

Pr(Y  <  y)  =  «[{(y/ic1)  °  -  V1  (hQ) } /a  (hQ)  ] ,  (4) 

2 

where  a  (h^)  *  ^2  C^0 )  Is  given  by  (2). 

It  is  well  known,  e.g.,  see  Anderson  (1958)  or  Srivastava  and 
Khatri  (1979),  that  for  many  likelihood  ratio  statistics  A  appear¬ 
ing  in  multivariate  analysis  under  the  normality  assumption, 

2  /N 

U  *  A  is  distributed  as  a  product  IIX^  of  independent  beta  var¬ 
iates  X±,  1*1,  2,  ...  k,  distributed  according  to  B(X±;  a^,  b^)  , 

where  N  Is  the  number  of  observations.  Equivalently,  we  have 
k 

-log  U  *  E  (-log  X.)  in  distribution.  Now,  It  can  be  shown  that, 
i-1 

as  a^^  and  b^-*30,  -log  X^  converges  in  law  to  normality.  Hence,  it 
is  possible  to  construct  a  normal  approximation  for  U  as  described 
above.  Towards  this  end  we  need  the  cumulants  of  -log  X^.  The 
moment  generating  function  of  -log  X^  is  easily  seen  to  be  M(t)  - 
B(a^-t,b^)/B(a^,b^) .  Hence,  the  cumulant  generating  function  is 
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K(t)  -  log{r(a±+b1)/r(ai)}  -  log{r(a1+b1-t)/r(ai-t)}. 
Differentiating  and  using  'ffZ)  *  log  T(Z),  the  cumulant  of 


-log  X.  is 


Crl  -  (-l)r{f(r’1)(a1)  -  'f(r_1)(ai+b1)}. 


But  ¥  (Z)  -  -  Z  (Z+j)  f  giving 
j-0 


Cri  *  C'-DW  Z  +  2  ((ai+j)“r  -  (ai+j+v)*r}  ]  ,  (5) 

j-0  j»m 

where  m  denotes  the  largest  integer  in  b^  and  v  *  b.^  -  m.  The 
cumulants  of  -log  Uf  obtained  using  (5/  are* 

k  m-1 

*  (U')  -  (r-1) ! [  Z  Z  (a.+j)‘r 
r  i-1  J-0  1 

k  00 

+  Z  Z  {(a^j)-*  -  (a1+j+v)"r}].  (6) 

i-1  j-m 

If  b^  is  an  integer  then  the  second  sum  in  (6)  vanishes  and 

k  m-1 

K  (O')  -  (r-1)’  I  Z  (a,+j )~r.  (7) 

r  i-1  j-0  1 

From  (7)  we  observe  that  as  either  k  or  b^  or  both  -*00,  diverges, 
but  <r>  r  >  1,  are  bounded.  That  is,  <J>r  -  /k^-*-0  .  Hence,  it  is 

possible  to  construct  the  normal  approximation  to  the  distribution 
of  A  as  described  above.  Thus,  from  (4)  we  get 

h  t 

Pr (A  >  X)  =  *[{(X,/k1)  °  -  V1(ho)}/O(ho)],  (8) 

where  X'  -  -  2(log  X)/N.  The  100(l-a)ttl  percentile  A^_a  can  be 
approximated  as 


Al-a  ^  +  W5 


(9) 


where  Z  denotes  the  100at^1  percentile  of  the  standard  normal 
ct 

variate • 

3.  TWO  APPLICATIONS  IN  MULTIVARIATE  ANALYSIS 
The  normal  approximation  derived  in  the  previous  section  is 
now  illustrated  and  later  examined  in  the  context  of  the  multi¬ 
variate  problems  of  testing  independence  between  two  sets  of  nor¬ 
mal  variates  and  testing  the  sphericity  hypothesis. 

3.1  Independence  Between  Two  Sets.  Let  Xf  *  (X^,  X2*  •••> 

X  )  and  Y 1  *  (Y- ,  Y9,  ...»  Y  ),  p-  <  p9,  p-  +  p9  *  p,  be  jointly 
p^  —  ±  Z  $2  ***  Z  X  z 

normally  distributed  with  Var(X)  *  E^,  Var(Y)  *  ^22  an<*  Cov(x»^) 

■  ^12*  The  hypothesis  of  Independence  between  X  and  Y  is 

■  0.  If  S  is  the  usual  estimate  of  £  based  on  a  sample  of  size  N, 
then  the  likelihood  ratio  statistic  for  testing  Hq  is 

a  -  [|s|/<|s11||s22|)]n/2, 

where  S^  and  S^  are  the  submatrices  of  S  corresponding  to  E^ 
and  ^22*  respectively.  The  exact  distribution  of  the  statistic  A 
is  given,  and  tabulated  for  some  values  of  p^  and  p^>  by  several 
authors  (e.g.,  see  Rrishnaiah  1979  and  Consul  1967a).  Among  var¬ 
ious  approximations  proposed  for  the  null  distribution  of  A  two 
are  well  known  and  widely  used  in  statistical  packages  such  as 
BMDP  (see  Engelman,  et  a_l.  ,  1977),  These  are  (i)  the  chi-square 
series  approximation  due  to  Box  (1949)  and  (ii)  the  F  approxi¬ 
mation  due  to  Rao  (1948). 

Box- Approximation.  Let  w  *  p^p2>  m  "  N  -  (p^+P2+3)/2, 

Y 2  *  w(p2+p2-5)/48,  y4  Y^2  +  w{3(p*+p2>  +  10w2  -  50(p2+p2) 


+  159} /1920 .  Then 


Pr(-m  log  0  z) 

"  1  z)  +  Y2(Pr(x^+4  <  z)  -  Pr(x2  1  * )>/» 2 

+  [Y4(Pr(x^.8  <  z)  -  Pr(x*  <  z)}  -  Y2{Pr(Xw+4  -  z) 

-  Pr(x^  <  z))]/mA  +  0(N'6),  (10) 

w  — 

2 

where  denotes  a  chi-square  variable  with  k  degrees  of  freedom 
and  U  -  A2/N. 

Rao- Approximation.  Let  m*  -  N  -  (pj+p2+3)/2,  L  *  (P1?2  “  2)/4, 
s  ■  (P^P2  ~  4)/(p2+p2  -  5)}.  Then, 

Q  -  (m’s  -  2L)(1  -  U1/S)/(PlP2U1/S) ,  (11) 

has  an  F-dis tr ibu tion  with  p^p£  and  m*s  “  2L  degrees  of  freedom* 

Now,  it  is  well  known  that  (e.g.,  see  Anderson,  1958,  p.  236) 

under  Hq  the  likelihood  ratio  statistic  A  satisfies  the  equivalence 
2  /N 

A  *  U  ■  HX±  in  law,  where  X^(I  *  1,  2,  .  ..,  P2)  are  indepen¬ 
dently  distributed  according  to  beta  distributions  B(X^;  (N  -  p^ 

-  i)/2,  p^/2).  The  normal  approximation  developed  in  the  previous 
section  can  be  specialized  in  this  case  by  taking  k  ■  " 

(N  -  p^  -  i)/2,  bi  ■  p^/2  in  the  expressions  (6)  for  the  cumulants, 
(8)  for  the  probabilities,  and  (9)  for  the  percentiles  of  A. 

3*2  Testing  the  Sphericity  Hypothesis.  Let  X^ ,  X^ »  ...»  XN 
be  a  random  sample  from  a  p-variate  normal  population  with  mean 
y  and  covariance  matrix  E.  The  hypothesis  that  the  p  components 

of  the  random  vector  X  are  independent  with  the  same  variance  i*e*, 

2  2 

Hq:E  -  a  1^,  a  >  0  unknown,  is  known  as  the  sphericity  hypothesis. 
The  hypothesis  also  arises  in  the  analysis  of  data  from  experiments 
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consisting  of  repeated  measurements.  In  these  experiments  ,  the 

measurements  on  a  subject  are  assumed  to  have  compound  symmetry, 

i.e.,  have  the  same  variances  and  same  correlations.  The  problem 

2 

of  testing  the  hypothesis  of  compound  symmetry  Hq:I  *  a  (pi  + 
(l-p)I)  for  the  covariance  structure  of  (p+1)  repeated  measure¬ 
ments  Y  can  be  reduced  to  the  sphericity  hypothesis  by  an  ortho¬ 
gonal  transformation  Y’-*Yf (l//(p+l) :T^)  where  1  is  the  vector  of 
l*s.  Y  satisfies  compound  symmetry  if  and  only  if  X  »  T^Y  satis¬ 
fies  the  sphericity  hypothesis.  The  likelihood  ratio  criterion 
for  the  sphericity  hypothesis  was  proposed  by  Mauchly  (1940)  as 

U  -  A2/N  -  I S I {(trS)/p}~p, 

where  S  is  the  covariance  matrix  of  the  sample  of  size  N.  He 
also  derived  its  null  distribution  for  p  =  2.  The  exact  null  dis¬ 
tribution  of  U  for  p  ~  3,  4  and  6  was  obtained  by  Consul  (1967b). 
The  51  and  1%  points  for  p  «  4(1)10  were  given  by  Nagarsanker  and 
Pillai  (1973).  The  series  approximation  due  to  Box  can  be 
expressed  in  this  case  as  follows. 

2 

Box- Approximation.  Let  e  *  p(p+l)/2  -  1,  f  -  n  -  (2p  +  p+2) 

/ ( 6p)  and  g  =  (p+2) (p-1) Cp-2) (2p3+6p2+3p+2)/(288p2)  for  n  -  N-l. 
Then, 

Pr(-  f  log  U  <_  z)  -  Pr(x2  <  z)  +  g{Pr(x^+4  <.  z) 

-  Pr(Xg  1  z))/f2  +  0(f-3).  (12) 

It  is  well  known  (e.g.  see  Sri vastava  and  Khatri,  1979)  that 

the  distribution  of  U  under  H  is  the  same  as  that  of  the  product 

o 

nxt,  where  X^  (i  ■  1,  2,  ...,  p-1)  are  independent  beta  random 
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variables  distributed  according  to  B(x^;  (n-i)/2,  i (p+2 ) / ( 2p) } , 

n  *  N  -  1*  Again,  we  can  obtain  the  cumulants  of  U*  ■ 

-  2 (log  U)/N  using  (6)  with  a±  -  (n-i)/2,  b±  -  i(p+2)/(2p)  and 
k  *  p  -  1.  Hence,  the  probabilities  and  the  percentiles  of  the 
likelihood  ratio  A  may  be  obtained  from  (8)  and  (9)f  respectively. 

4.  NUMERICAL  COMPARISONS 

The  quality  of  the  normal  approximations  for  the  two  multi¬ 
variate  likelihood  ratio  statistics  discussed  in  the  previous 
section  and  the  other  two  approximations,  examined  by  comput¬ 

ing  the  probabilities  corresponding  to  the  tabulated  percentiles 
of  the  statistics.  Thus,  in  the  case  of  the  null  distribution  of 
A  for  testing  independence,  the  approximations  due  to  Box  (10), 
due  to  Rao  (11),  and  the  normal  approximation  given  in  Section 
(3.1)  were  used  to  compute  the  probabilities  corresponding  to  all 
5 Z  and  1%  points  of  A  given  in  Pearson  and  Hartley  (1972,  p.  99 
and  333).  Similarly,  in  case  of  the  sphericity  problem,  all  per¬ 
centiles  given  by  Nagarsanker  and  Pillai  (1973)  were  used  to 
examine  the  approximation  due  to  Box  given  by  (12)  and  the  relevant 
normal  approximation.  In  both  cases,  the  series  approximation  due 
to  Box  was  used  in  two  steps:  1)  only  the  first  term;  and  2)  all 
terms  given  In  (10)  and  (12).  Also  the  percentiles  approximated 
using  the  normal  approximations  were  compared  with  the  competing 
approximations  using  the  first  term  of  the  Box  series  and  the  F- 
approximat ion.  A  selection  of  errors,  i.e.,  (Approximation  - 
Exact  value)X10^,  in  various  cases  is  presented  in  Tables  1  and  2. 
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Conclualons .  Let  New,  Rao,  Box  1  and  Box  3  denote  the  nor¬ 
mal  approximation,  the  F-approximat ion  due  to  Rao,  the  first  term 
approximation  due  to  Box  and  the  three  term  approximation  due  to 
Box,  respectively.  From  Tables  1  and  2  it  may  be  observed  that 
(i)  Rao,  Box  1  and  Box  3  have  errors  in  second  through  fifth 
decimal  place,  they  are  especially  large  for  small  N  and  decreas¬ 
ing  rapidly  as  N  increases.  The  normal  approximation  has  errors 
in  fourth  or  fifth  decimal  place.  (ii)  As  p^,  p^  or  p  increases, 
errors  due  to  Rao,  Box  1  and  Box  3  increase  while  those  due  to  the 
normal  approximation  either  decrease  or  maintain  the  same  level. 
Overall,  the  normal  approximation  is  superior  for  small  N  and  is 
comparable  with  the  others  when  the  N  is  large. 
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TABLE  1.  Errors  of  the  Approximations  for  the  Likelihood  Ratio  Statistic 

for  Testing  Independence  Between  Two  Sets 

a  -  .05 

PI  P2  N  X  ERRORS* 

PERCENTILES  PROBABILITIES 


NEW 

RAO 

B0X1 

NEW 

RAO 

BOX1 

BOX3 

8 

12 

0.00001 

0 

0 

67 

-19 

-1537 

-4991 

-4744 

8 

19 

0.04107 

4 

7 

698 

-17 

-24 

-1776 

-50 

22 

30 

0.00107 

0 

1 

356 

0 

-143 

-4936 

-3893 

22 

37 

0.01620 

0 

0 

760 

-4 

0 

-3842 

-725 

8 

18 

0.00217 

0 

5 

161 

2 

-218 

-3340 

-485 

8 

25 

0.03411 

1 

3 

352 

-8 

-17 

-1355 

-16 

16 

26 

0.00019 

0 

0 

50 

12 

-420 

-4759 

-2732 

16 

33 

0.00568 

0 

1 

213 

-6 

-48 

-3160 

-316 

8 

19 

0.00021 

0 

1 

37 

-9 

-618 

-4225 

-1437 

8 

23 

0.00356 

0 

3 

133 

-3 

-122 

-2657 

-199 

10 

21 

0.00007 

0 

0 

18 

25 

-746 

-4531 

-2059 

10 

25 

0.00157 

0 

2 

81 

9 

-155 

-3153 

-358 

a  - 

.01 

8 

12 

0.00000 

0 

0 

17 

6 

-478 

-999 

-997 

8 

19 

0.02261 

0 

7 

514 

0 

-9 

-487 

-29 

22 

30 

0.00043 

0 

0 

209 

0 

-44 

-997 

-935 

22 

37 

0.00990 

-2 

-1 

572 

6 

5 

-872 

-303 

8 

18 

0.00085 

0 

3 

85 

3 

-72 

-803 

-223 

8 

25 

0.02086 

0 

6 

270 

-2 

-11 

-374 

-15 

16 

26 

0.00007 

0 

0 

25 

-10 

-143 

-985 

-784 

16 

33 

0.00327 

0 

0 

147 

3 

-11 

-753 

-147 

8 

19 

0.00007 

0 

0 

16 

1 

-192 

-932 

-513 

8 

23 

0.00174 

0 

2 

81 

-1 

-42 

-665 

-102 

10 

21 

0.00002 

0 

0 

8 

3 

-234 

-965 

-659 

10 

25 

0.00075 

0 

1 

48 

5 

-48 

-756 

-164 

TABLE  2.  Errors  of  Approximations  for  the  Likelihood  Ratio  Statistic  for 

Testing  Sphericity 

a  -  .05  a  -  .01 

P  N  X  ERRORS*  ERRORS* 


PERC. 

PROB. 

X 

PERC 

J  . 

PROB. 

NEW 

B0X1 

NEW 

B0X1 

B0X2 

NEW  B0X1 

NEW 

BOXl 

B0X2 

4 

10 

0.09739 

16 

467 

-21 

-539 

-57 

0.05010 

-22 

382 

10 

-173 

-25 

15 

0.25350 

35 

236 

-29 

-190 

-10 

0.17210 

-41 

267 

9 

-62 

-4 

20 

0.37720 

43 

120 

-33 

-97 

-5 

0.28670 

-44 

166 

8 

-31 

-1 

30 

0.53900 

46 

37 

-38 

-40 

-2 

0.45310 

-38 

76 

7 

-12 

0 

5 

10 

0.03110 

1 

444 

-7 

-1171 

-236 

0.01361 

-5 

281 

7 

-350 

-96 

15 

0.13780 

7 

334 

-10 

-393 

-29 

0.08685 

-19 

287 

7 

-124 

-14 

20 

0.24820 

17 

201 

-17 

-194 

-7 

0.17970 

-23 

188 

6 

-63 

-5 

30 

0.41640 

28 

79 

-27 

-78 

-2 

0.34020 

-29 

58 

6 

-23 

0 

7 

10 

0.00094 

0 

104 

20 

-3352 

-1789 

0.00025 

0 

43 

2 

-822 

-566 

15 

0.02712 

0 

297 

-4 

-1158 

-228 

0.01444 

-3 

207 

4 

-330 

-86 

20 

0.08446 

5 

282 

-12 

-562 

-63 

0.05514 

-5 

234 

3 

-165 

-24 

30 

0.21780 

15 

163 

-21 

-216 

-13 

0.16770 

-4 

153 

1 

-65 

-5 

♦Error  in  prob.  ■  (Approx,  value  -  a)xl0^  and  in  perc.  (Approx,  value  -  X)xl0^. 


-12- 


REFERENCES 

Anderson,  T.  W.  (1958).  An  Introduction  to  Multivariate  Statis¬ 
tical  Analysis.  John  Wiley  and  Sons,  New  York. 

Bahadur,  R.  R.  (1967).  An  optimal  property  of  the  likelihood 
ratio  statistics.  Proceedings  of  Fifth  Berkley  Symposium 
Mathematical  Statistics  and  Probability.  Univ.  of  California 
Press,  13-26. 

Bartlett,  M.  S.  (1937).  Properties  of  sufficiency  and  statistical 
tests.  Proceedings  of  Royal  Statistical  Society,  A,  160,  268. 

Bishop,  D.  J.  (1939).  On  the  comprehensive  test  of  the  homogenec- 
ity  of  variances  and  covariances  in  multivariate  problems. 
Biometrika  a  31,  31-55. 

Box,  G.  E.  P.  (1949).  A  general  distribution  theory  for  a  class 
of  likelihood  criteria.  Biometrika  t  36,  317-46. 

Consul,  P.  C.  (1967a).  On  the  exact  distribution  of  likelihood 
ratio  criteria  for  testing  independence  of  sets  of  variates 
under  the  null  hypothesis.  Annals  of  Mathematical  Statistics, 
38,  1160-9. 

Consul,  P.  C.  (1967b).  On  the  exact  distribution  of  the  criterion 
W  for  testing  sphericity  in  a  p-variate  normal  distribution. 
Annals  of  Mathematical  Statistics.  38,  1170-4. 

Engelman,  L.,  Frane,  J.  W. ,  and  Jennrich,  R.  I.  (1977).  Biomedical 
Computer  Programs  P-Series.  Univ.  of  California  Press, 

Berkley,  California. 

Hsieh,  H.  K.  (1979).  On  asymptotic  optimality  of  likelihood  ratio 
tests  for  multivariate  normal  distributions.  Annals  of 
Statistics ,  7,  592-8. 

Krishnaiah,  P.  R.  (1979).  Some  Recent  Developments  on  The  Real 
Multivariate  Distributions.  North  Holland  Co.,  New  York. 

Mann,  H.  B.  and  Wald,  A.  (1943).  On  stochestic  limit  and  order 
relationships.  Annals  of  Mathematical  Statistics,  14,  217-26. 

Mauchly,  J.  W.  (1940).  Significance  test  for  sphericity  of  a  nor¬ 
mal  n-variate  distribution.  Annals  of  Mathematical  Statistics. 
11,  204-9. 

Nagarsanker,  B.  N.  and  Pillai,  K.  C.  S.  (1973).  Distribution  of 
the  sphericity  test  criterion.  Journal  of  Multivariate 
Analysis ,  3,  226-35. 


-13- 


Nayer,  P.  P.  N.  (1936).  Statistical  Research  Memorandum,  1,  38. 

Neyman,  J.  and  Pearson,  E.  S.  (1931).  Bulletin  International 
Academy  Cracovie,  A,  460. 

Pearson,  E.  S.  and  Hartley,  H.  0.  (1972).  Biometrika  Tables  for 

Statisticians ,  2.  Cambridge  Univ.  Press,  New  York 

Rao,  C.  R.  (1948).  Tests  of  significance  in  multivariate  analysis. 
Biometrika ,  35,  58-79. 

Srivastava,  M.  S.  and  Khatri,  C.  G.  (1979).  An  Introduction  to 
Multivariate  Statistics,  North  Holland  Co.,  New  York 


SECURITY  CLASSIFICATION  OF  THS  PAGE  fHbrn  fUna  hnte 


REPORT  DOCUMENTATION  PAGE 


RKAI)  IN S'l  ,UT‘  l  ION!. 

UM  OkK  roMi'i.K  i  i.nt;  r>km 


RE  PORT  NUMBER 


AFOSR-TR-  8  1-0180 

TITLE  fond  Subtitle) 

A  NORMAL  APPROXIMATION  FOR  1HE  MULTIVARIATE 
LIKELIHOOD  RATIO  STATISTICS 


AUTMORfs.) 

Govinci  S.  Muclholkar  and  Madhusudan  C.  Trivedi 


9  performing  organization  name  and  address 

University  of  Rochester 
Department  of  Statsitics 
Rochester,  N.  Y.  14627 


11.  controlling  office  name  anO  address 

Air  Force  Office  of  Scientific  Research/NM 
Bolling  AFB,  Washington,  DC  20332 


2  GOVT  ACC  CSSION  NO  I  3  RF  CIP1  f.N  T’S  •'  A  T  A  t  OG  NUM11 1  R 


S  TYPE  OF  REPORT  A  P  L  RIOD  COVE  R  F  I 


Interim 


6  PERFORMING  o^g.  rlport  NUMBER 


8  CONTRACT  OR  GRANT  NUMBER'  s' 


AFOSR- 77-3360 


10  program  element,  project,  TA,K 
AREA  6  WORK  UNIT  NUMBERS 


61102F  2304/A5 


12.  REPORT  DATE 


13  NUMBER  OF  PAGES 


MONITORING  AGENCY  NAME  ft  AODRESSfi/  different  from  Controlling  Office)  15.  SECURITY  CLASS,  (of  this  report) 

UNCLASSIFIED 

D £ C L  A SS»  fT< “AT |ON  OO  W N  GR  A O l  N 
SCH  EDULE 


DISTRIBUTION  STATEMENT  (of  this  Report) 


Approved  for  public  release;  distribution  unlimited 


DISTRIBUTION  STATEMENT  (of  the  abstract  entered  m  Rlnck  20.  it  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


19  KEY  WORDS  (Continue  on  reverse  side  it  necess/ irv  nnd  identitv  by  block  number) 


BSTRACT  (Cmiflmie  nfi  rpvrfsc  aide  If  necessary  and  Identify  bv  block  number) 

{r0R  MANY  MULTIVARIATE  HYPOTHESES,  UNDER  THE  NORMALITY  ASSUMPTIONS,  THE 
LIKELIHOOD  RATIO  TESTS  ARE  OPTIMAL  IN  THE  SENSE  OF  HAVING  MAXIMAL  EXACT  SLOPES. 
THE  EXACT  DISTRIBUTIONS  NEEDED  FOR  IMPLEMENTING  THESE  TESTS  ARE  COMPLEX  AND 
THEIR  TABULAIT0N  IS  LIMITED  IN  SCOPE  AND  ACCESSIBILITY.  IN  THIS  PAPER,  A 
METHOD  OF  CONSTRUCTING  NORMAL  APPROXIMATIONS  TO  THESE  DISTRIBUTIONS  IS 
DESCRIBED,  AND  ILLUSTRATED  USING  THE  PROBLEMS  OF  TESTING  SPHERICITY  AND 
INDEPENDENCE  BETWEEN  TOO  SETS  OF  VARIATES.  THE  NORMAL  APPROXIMATIONS  ARE 
COMPARED  WITH  WELL  KNOWN  COMPETING  APPROXIMATIONS  AND  SEEN  TO  FARE  WELL 


*-*?“**V 


V  - 


